Temporal self-restoration of compressed optical filaments 
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We numerically investigate the propagation of a self- compressed optical filament through a gas- 
glass-gas interface. Few-cycle light pulses survive a sudden and short order-of-magnitude increase 
of nonlinearity and dispersion, even when all conservative estimates predict temporal spreading or 
spatial breakup. Spatio-temporal distortions are shown to self-heal upon further propagation when 
the pulse refocuses in the second gas. This self-healing mechanism has important implications for 
pulse compression techniques handled by filament at ion and explains the robustness of such sources. 



The physics of ultrashort optical filaments has re- 
cently received increased attention as it enables self- 
compression of intense laser pulses into the few-cycle 
regime [l], 0, uj m B B 0- Femtosecond pulses with 
input powers close to the self-focusing threshold, P cr , in- 
deed produce a single filament in which the spectrum 
broadens and the pulse duration decreases upon propa- 
gation. Experimental 0, 0] as well as theoretical work 
@, 0] indicated that self-compression occurs in a limited 
parameter range, and small deviations from these ideal 
conditions may either cause spatial breakup into multi- 
ple filaments or temporal splitting of the pulse. For ev- 
idencing self-compression, experiments currently employ 
Brewster-windowed cells filled with noble gas at pres- 
sures controlling the ratio of input power over critical. 
Diagnostics are furthermore placed several tens of cm be- 
yond the sample, making the pulse undergo a last stage 
of atmospheric propagation. A crucial point is then to 
preserve the main characteristics of the compressed pulse 
when it exits the output window. Several points concern- 
ing this issue, however, remain unclear. First, the pulse 
entering the glass sample undergoes a nonlinear index be- 
ing three decades higher than in the gas. This renders the 
robustness of the exiting optical structure questionable as 
a light beam with several thousands of critical powers is 
expected to immediately break up by modulational insta- 
bility @,!M[ni. 

Second, the much stronger dispersion 
of glass should severely alter the temporal compression. 
For instance, a 1 mm thick sample with group- velocity 
dispersion (GVD) coefficient ~ 370 fs 2 /cm at 800 
nm can double the duration of a 8 fs transform-limited 
Gaussian pulse [l^ . 

In the following we will show that despite these appar- 
ent sources of instability, filamentary propagation holds 
a hidden self-healing mechanism for ultrashort pulses. 
Comparable phenomena have already been observed in 
several physical areas, e.g., in solid-state dynamics on 
atomic scales [l3| and in plasma physics under mag- 
netic confinement P3], explaining stabilization mecha- 
nisms in apparently unstable systems. In optics, spatial 
self-restoration of Bessel beam profiles [l5| and of fila- 
ment cores [H, have also been reported. In contrast 
to all these previous observations, we report on simul- 
taneous spatio-temporal self-restoration of light bullets 



having experienced a perturbation that should have im- 
mediately caused spatial and temporal decay of the fil- 
ament. In this Letter we identify the changes in a self- 
compressed filament when it leaves a cell filled with ar- 
gon, propagates in a mm long silica window, and finally 
travels through the atmosphere. We show that glass- 
dispersion acts jointly with the Kerr response to widely 
stretch the pulse in time. Unexpectedly, the filament 
eventually recompresses in air to restore temporal profiles 
with few-cycle durations. This self-healing mechanism 
has not an anecdotic character, as it can be responsible 
for the remarkably short pulse durations seen in many 
self-compression experiments. 

Propagation equations model the forward component 
of the pulse envelope, S(x 1 y j _zA) 1 coupled with the free 
electron density, p(x, y, z, t) [l8[ 
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where / = |£| 2 , z is the propagation variable and t the 
retarded time in a reference frame moving with group ve- 
locity 1/k^ at laser mid-frequency coo; k(u) = n(uj)uj/c, 
ko = k(uo), no = n(ouo) and T = 1 + (i/uo)d t [19]. V = 
^Z n>2 (^ n - > /n\)(idt) n is the dispersion operator formally 
involving the derivatives k^ = d n k/dcj n \ UJo . Linear dis- 
persion curves n(u) are taken from Refs. [l2|, [20|, Hl| 
for argon, silica and air, respectively. At 800 nm, no 
takes the value 1 for gases and 1.45 for glass. The 
leading GVD term in V [k^] is about three orders of 
magnitude larger in glass than in the two gases; p c = 
1.73 x 10 21 cm -3 is the critical plasma density. The 
nonlinear Kerr response contains only an instantaneous 
contribution, TZ(t) = when the pulse self- focuses 

in argon. It includes a Raman-delayed contribution, 



K(t) = (1 - f)S(t) 
ratio / in silica [22I , 



|2J and in air 



-^ TR sm(u R t), with 
If. Plasma sources 



are driven by the PPT ionization rate W(I) [25[ in argon 
and air, and by the Keldysh rate for crystals in glass [261 ]. 
In the argon cell we assume a uniform pressure of 0.5 bar, 
leading to a nonlinear Kerr index n 2 = 5 x 10 -20 cm 2 /W 
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at 800 nm. The Kerr index for glass is n 2 = 3.2 x 10 -16 
cm 2 /W and = 2.5 x 10 -19 cm 2 /W in air at atmo- 
spheric pressure . All other parameter values used in 
this paper, including ionization potentials avalanche 
cross-sections <r, neutral densities p n t, and electron re- 
combination time for silica r rec , can be found in Ref. Q 
(argon) and Ref. [18[ (glass, air). Simulations have been 
performed in full 3D (Vj_ = d 2 + d 2 ) and radially- 
symmetric (V^_ = r~ 1 d r rd r , r = \J x 2 + y 2 ) geometries, 
yielding analogous results. The starting pulse, being 
Gaussian, is focused by a lens of focal length / = 50 cm 
into an argon cell with a 1.5 m maximum length. Input 
waist at 1/e 2 and full width at half-maximum (FWHM) 
duration are wq = 500 /im and To — 30 fs. The input 
power Pi n is equal to one critical power in argon, i.e., 
P cr Ao/27rnon2 = 20.4 GW. The following analysis is 
divided into three steps devoted to pulse propagation in 
argon, inside the silica window, and propagation in air. 

1. Filamentation in argon. With one critical power, 
the pulse develops a light bullet dynamics 0], charac- 
terized by a primary focusing sequence, followed by a 
second, limited refocusing [Fig. DJa)]. Minimum FWHM 
duration is attained when the plasma relaxes, and the 
filament is compressed to about 6 fs at z ~ 0.7 m [Fig. 
GJb)]. The three subplots of Fig. [IJc) show intensity 
distributions at distances z = 0.5, 1 and 1.5 m where 
the on- axis fluence {T = J \£ \ 2 dt) reaches 1.7 J/cm 2 , 0.2 
J/cm 2 , and 0.05 J/cm 2 , respectively. These distances de- 
fine the effective length of the argon cell and the position 
of the exit glass window. 

2. Crossing the silica window. At high enough flu- 
ences > 0.1 J/cm 2 , only a fraction of the forward pulse 
can be transmitted through a silica surface. According to 
[27[ | , only 20 % of the forward pulse is expected to trans- 
mit the glass surface at z = 0.5 m, whereas this fraction 
increases to 90 % and 100 % for z = 1 and 1.5 m, respec- 
tively. Conditions beyond z = lm(/<2x 10 13 W/cm 2 ) 
therefore avoid damage of silica [Hf. Figure [2](a) summa- 
rizes our results. Output profiles of Fig. [TJc) multiplied 
by the previous transmission rates are used as initial con- 
ditions of Eqs. ([I])-© solved for silica. 3D and radially- 
symmetric simulation results almost superimpose; hence, 
no severe spatial distortions affect the intensity profiles 
in the (x,y) plane. The pulse captured at z = 0.5 m 
(pulse I) begins to diffract over 1 mm, but it refocuses 
and collapses just afterwards. Although experimentally 
irrelevant, this extreme configuration already points out 
the importance of the first few mm of propagation dom- 
inated by pulse dispersion. Triggered by beam fluences 
> 1.5 J/cm 2 in argon, the surface plasma created by 
pulse I should, however, act as a filter transmitting only 
the leading edge of the pulse [29[, thereby rendering the 
underlying assumption of a fixed temporal profile invalid. 
For pulses II and III having much lower fluences, assump- 
tion of a fixed profile still holds. In the following, we thus 
concentrate on these last two configurations only. 

When the window is positioned at z = 1 m (pulse II), 
again dispersion first prevails, but does not prevent beam 
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FIG. 1: (a) Peak intensity (measured on left axis) and elec- 
tron density (resp. right axis) for a 500 /xm-waisted, 30 fs 
Gaussian pulse focused with / = 50 cm inside a 1.5 m long 
cell of argon computed in 3D geometry (solid curves) and in 
radial symmetry (dotted curves). 3D simulations involve a 
5% amplitude random noise in the input, (b) Temporal on- 
axis dynamics, (c) Intensity profiles in the (x,t) plane at the 
three positions selected for the exit glass window. 
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FIG. 2: (Top) Peak intensity of the filament exiting the argon 
cell at (I) z = 0.5 m, (II) z = 1 m and (III) z = 1.5 m, and 
propagating into the silica window over mm distances. (Mid- 
dle and bottom) Rows (II) and (III) detail the (x,t) evolution 
of the pulse profiles II and III inside the exit window. 



collapse over 5 mm long paths. The input beam has a 
waist of 240 jiva and conveys time slices with a maxi- 
mum power of - 7100 P cr , where Pg = 2.19 MW denotes 
the critical power in glass. In contrast, the low- intensity 
pulse III disperses rather slowly over the same distance 
and beyond. The maximum input power, however, is still 
about 8000 P cr , as the beam waist entering the sample 
has been enlarged to ~ 0.6 mm. Rows (II) and (III) of 
Fig. [2] depict the evolution of pulses II and III in the 
(x, t) plane. Pulse II enters the dielectric with time ex- 
tent of ~ 8 fs. After 0.5 mm in silica, this pulse decays 
into a broader structure of ~ 33 fs, forming a pedestal 
that extends deep into the trailing region. After 1 mm, 
its duration increases to 50 fs. Pulse III first exhibits a 
duration of 13 fs, which then broadens to 28 fs over 0.5 
mm in glass and more afterwards. 

It is important to note the complete absence of mul- 
tifilamentation (MF) in our simulations. Estimates from 
modulational instability theory Q, however, predict an 
onset for MF and subsequent collapse upon short dis- 
tances, A^mf — ^o^cr/Ao^max ~ 0.1 mm, for peak in- 
tensities Io > 3 TW/cm 2 impinging on the glass sur- 
face. In fact, as GVD rapidly damps the peak intensity, 
the amplification of periodic spatial modulations lead- 
ing to MF is practically prevented. Here the plasma 
density plays a negligible role in the nonlinear index, as 
An N L - n 2 /max-Pmax/2n pc ~ l-5x 10~ 3 . We also nu- 
merically checked that high-order dispersion, pulse steep- 
ening and Raman scattering are of minor influence. 

Besides, it turns out that both GVD and the Kerr 
response are responsible for the strong increase of the 
pulse duration that exceeds predictions by the well- 
known linear formula for Gaussian pulses, T\\ n (z) 
rov 7 ! + 16(ln2) 2 (/c( 2 )Az/r 2 ) 2 0. Denoting the spa- 
tial and temporal radii of Gaussian pulses by R(z) and 
T(z) = r(z)/2Vln 2, respectively, a two-scaled varia- 
tional principle applied to the reduced equation id z i/j + 
Vj_-0 — Sofa + l^l 2 ^ = can ne fp us understand these 
peculiar dynamics through the system [3(j 

±&R„ = 1-JL ; l I * T „ = 6(6+^ i ). (3) 

In Eq. ©, p = Pin/ Per and S = 27m i? 2 il /c (2) /XoT^ is the 
normalized GVD coefficient for the input filament with 
size R n \ and temporal radius T n \. 5 takes very large val- 
ues > 200 when using durations and waists of pulses II 
and III. Since the growth of T(z) is driven by S 2 (dis- 
persion) and S x p (Kerr response), the beam collapse 
characterized by the vanishing R(z) — > at finite z can 
efficiently be delayed in the limits J, p ^> 1. Indeed, the 
more T(z) increases, the less R(z) decreases. Moreover, 
Eq. ([3]) illustrates that both GVD and cubic nonlinear- 
ity contribute to enlarge the pulse duration. Typically, 
FWHM durations of ~ 25 fs are attained over 0.5 mm by 
the initially 8 fs long pulse II, whereas linear dispersion 
alone allows a temporal broadening limited to ~ 10 fs. 
For this pulse, broadening is more than doubled by the 
Kerr response. Hence, GVD disperses the pulse while at 
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FIG. 3: On- axis temporal evolution of (a) pulse II and (b) 
pulse III in the atmosphere, after exiting a 0.5 mm thick win- 
dow, (c) Temporal profiles compressed in air at the distances 
z — 0.25 m (solid line) and 0.5 m (dash-dotted line) for pulse 
II. (d) Same for pulse III at z — 0.5 m (solid line) and z — 1 
m (dashed line). 



high powers the Kerr term stretches it in time. Combin- 
ing both effects causes significant pulse broadening. 

3. Air propagation. Because diagnostics are usually 
positioned tens of cm away from the cell, the beam finally 
propagates in atmosphere. Figures E9(a,b) show the tem- 
poral evolution of pulses II and III as they exit a 0.5 mm 
thick window and are fully transmitted in air. These two 
configurations are characteristic of filaments diagnosed 
either at moderately high intensities (II) or close to the 
end of the diffraction stage (III). Because these pulses 
contain time slices with maximum powers of ~ 3.5 P cr in 
air, they eventually refocus and thus recompress, which 
is the core of the self-healing process. Pulse II attains 
intensities less than 50 TW/cm 2 and triggers a minor 
plasma response. Pulse III involves much weaker inten- 
sities. The key-player of this process is the self-focusing 
mechanism alone, which amplifies the time slice with the 
largest power to the detriment of the others [18[. Col- 
lapse is not arrested by plasma generation, but instead by 
Raman-scattering. The Raman-delayed response indeed 
doubles the effective power for self-focusing and delocal- 
izes the pulse in time. At powers close to critical, this 
is sufficient to arrest the beam collapse. Figures [3^c,d) 
detail temporal profiles with FWHM durations of 10.3 fs 
at z = 0.25 m and 10.8 fs at z = 0.5 m for pulse II, and of 
16.6 fs and 18.6 fs at 0.5 m and 1 m for pulse III. These 
durations are relatively stable since they hold over sev- 
eral tens of cm. They remain in the order of magnitude 
of those reported in the experiments Q. 

Figure [4] summarizes the history of on-axis spectra 
around the pump frequency (coo = 2.35 PHz) along the 
propagation steps 1 to 3. Figure Ufa) describes pulse 
II, passing through the glass window and transmitted to 
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FIG. 4: On-axis spectra for (a) pulse II and (b) pulse III. Red 
curves refer to the dynamics after compression in argon. Blue 
curves show that in the window at z = 0.5 mm. Black curves 
refer to air propagation with data collected at z = 0.5 m. 

air. The first focusing event in argon creates red- and 
blue-shifted spectral wings caused by self-phase modula- 
tion (red curve). Next, GVD drives the propagation in 
the dielectric, which lowers the saturation intensities and 
weakens self-phase modulation, such that the blueshifted 
wing drops out (blue curve). Figure 01(b) provides the 
same pieces of information for pulse III. At weaker in- 
tensities, the spectral evolution follows that of pulse II, 
although the properties expected from linear dispersion 
become more important and plasma turns off. Hence, 



on-axis spectra remain closer to those exiting argon. 

It is worth emphasizing that other compressor schemes 
based, e.g., on an Ar-glass-Ar interface, were also numer- 
ically simulated (not shown). Increasing the pressure to 
1 bar in the second argon cell guarantees a power above 
critical for pulse II, which then refocuses and follows a 
temporal evolution similar to that of Fig. [3ja). This 
underlines the generality of the self-healing mechanism. 

In conclusion, numerical simulations cleared up the dy- 
namics of self-compressed filaments when they cross a 
gas-glass-gas interface. After they are produced in a no- 
ble gas, they propagate over a few mm inside a silica 
window and run tens of cm in air before being measured. 
In the window, the interplay between dispersion and Kerr 
self-focusing inhibits multifilamentation and significantly 
increases pulse durations. Despite this large stretching 
in time, durations in the 10 fs range restore themselves 
upon further nonlinear propagation in air. This new self- 
healing mechanism acts both in space and time. Unlike 
previous scenarios where the perturbation is caused by 
loss mechanisms, such as raindrops, apertures or wire 
grids, the self-restoration is here induced by a discontinu- 
ity of the nonlinear properties inherent to the materials. 
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